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Using simulations of helically driven turbulence, it is shown that the ratio of kinetic to magnetic energy dissipation scales 
with the magnetic Prandtl number in power law fashion with an exponent of approximately 0.6. Over six orders of mag- 
nitude in the magnetic Prandtl number the magnetic field is found to be sustained by large-scale dynamo action of alpha- 
squared type. This work extends a similar finding for small magnetic Prandtl numbers to the regime of large magnetic 
Prandtl numbers. At large magnetic Prandtl numbers, most of the energy is dissipated viscously, lowering thus the amount 
of magnetic energy dissipation, which means that simulations can be performed at magnetic Reynolds numbers that are 
large compared to the usual limits imposed by a given resolution. This is analogous to an earlier finding that at small 
magnetic Prandtl numbers, most of the energy is dissipated resistively, lowering the amount of kinetic energy dissipation, 
so simulations can then be performed at much larger fluid Reynolds numbers than otherwise. The decrease in magnetic 

energy dissipation at large magnetic Prandtl numbers is discussed in the context of underluminous accretion found in some 
quasars. 
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1 Introduction 

The magnetic fields in astrophysical bodies often have a 
pronounced large-scale component that is associated with 
large-scale dynamo action. Examples are the cyclic mag- 
netic fields in late-type stars such as the Sun and the mag- 
netic spirals in many galaxies, including even irregular 
galaxies; see Beck et al. (1996) for a review. In addition, 
all observed magnetic fields also have a significant small- 
scale component that may either be the result of turbulent 
motions distorting the large-scale field, or, alternatively, it 
could be the result of what is known as small-scale dynamo 
action (Cattaneo 1999). 

Much of our knowledge about large-scale and small- 
scale dynamos has come from numerical simulations; see 
Brandenburg & Subramanian (2005) for a review. It is clear 
that, in order for simulations to approach an astrophysically 
interesting regime, one wants to make both the magnetic 
diffusivity and the kinematic viscosity as small as possible. 
This means that the magnetic and fluid Reynolds numbers 
should be as large as possible for a given numerical reso- 
lution, N 3 . The relevant criterion for sufficient numerical 
resolution is that the kinetic and magnetic energy spectra 
should develop an exponentially decaying dissipative sub- 
range at a wavenumber that is at least a factor of 10 be- 
low the Nyquist frequency, k-^ y = irN/L. In practice, for 
example, with a simulation at a resolution of 512 3 mesh 
points, one can hardly exceed values of the magnetic and 
fluid Reynolds number of about 500-700 (e.g., Branden- 
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burg 2009). However, as will be discussed in more detail 
in this paper, this empirical constraint on the resolution re- 
ally only applies if the ratio of magnetic and fluid Reynolds 
numbers is about unity. This ratio is also referred to as the 
magnetic Prandtl number, Pim, and there is hardly any sys- 
tem where this number is unity. In galaxies and galaxy clus- 
ters this number tends to be very large, while in stars and 
stellar accretion discs it is quite small. Also liquid metals 
used in laboratory experiments have small Pim- Therefore, 
much of what has been learnt from numerical simulations 
at Pim ~ 1 has to be re-examined in cases of low and high 
values of Pim- 

The purpose of this paper is to focus on the relative 
importance of viscous and ohmic dissipation rates at dif- 
ferent values of Pr^ . Often, viscous and ohmic dissipation 
are only treated "numerically" by making sure the code is 
stable. In such cases, viscosity and magnetic diffusivity are 
usually not even stated explicitly in the equations, suggest- 
ing that these terms are negligible and not important. This 
is of course not the case, as can be illustrated by consid- 
ering the case of quasars that belong to the most luminous 
objects in the sky. The discovery of the first quasar, 3C 273, 
is nicely explained by Rhodes (1978) in a popular maga- 
zine. Indeed, 3C 273, has about 2 x 10 12 times the lumi- 
nosity of the Sun and is indeed the brightest one in the sky. 
This quasar would not shine at all if it was not for the ef- 
fect of microphysical viscosity that leads to viscous dissi- 
pation. But how important is viscous dissipation compared 
with ohmic dissipation? In order to address this problem we 
need to understand the effects of both viscosity and mag- 
netic diffusivity in a turbulent system where the magnetic 
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field is self-sustained by dynamo action. In this paper we 
review briefly some recent work on dynamos in the regime 
of small PrM and turn then to the investigation of large PrM- 

2 Small magnetic Prandtl number dynamos 

In the last 6 years the issue of low magnetic Prandtl num- 
bers, PrM = v/rj, has become a frequently discussed topic 
in the dynamo community. This is the regime where the 
magnetic diffusivity r\ is large compared with the kinematic 
viscosity v. Already over a decade ago, Rogachevskii & 
Kleeorin (1997) noticed that for small-scale dynamos the 
critical value of the magnetic Reynolds number, ReM, for 
the onset of dynamo action should rise from a value around 
35 at PrM = 1 to values around 400 for small values of 
PrM- Here, ReM = u Ilns /r]kf is defined with respect to 
the wavenumber fcf of the energy-carrying eddies and the 
rms velocity, u rms . However, the result of Rogachevskii & 
Kleeorin was not widely recognized at the time. In 2004, 
simulation began to address this point systematically. Sim- 
ulations of Schekochihin et al. (2004) and Haugen et al. 
(2004) provided clear indications that Re M '* rises, and the 
results of Schekochihin et al. (2005) might have even sug- 
gested that the critical value of ReM for small-scale dynamo 
action might have become infinite for PrM ~ 0.1. 

Meanwhile, Boldyrev & Cattaneo (2004) provided an 
attractive framework for understanding this behavior. Given 
that the energy spectrum of the small-scale dynamo peaks 
at the resistive scale, which is the smallest possible scale at 
which the motions can still overcome resistive damping, one 
must ask what are the properties of the flow at this scale. 

In the original scenario of Kazantsev (1968), the small- 
scale dynamo works through a velocity field that is random, 
but essentially laminar and of large scale. In a simulation 
this can be realized by choosing a large magnetic Prandtl 
number, so the magnetic Reynolds number is much larger 
than the fluid Reynolds number. However, subsequent stud- 
ies show that small-scale dynamo action can also occur for 
PrM of order unity. Both for PrM = 1 and for PrM S> 1 
one finds that the spectral magnetic energy increases with 
wavenumber proportional to fc 3 / 2 . 

A qualitatively new feature emerges when PrM is small. 
In that case the wavenumber corresponding to the resistive 
scale decreases and lies in the inertial range of the turbu- 
lence. This property is crucial because in the inertial range 
the velocity field is "rough", i.e. over a spatial interval Sx 
the velocity difference Su = u(x + Sx) — u(x) scales like 
Su <~ Sx^ where £ < 1. Thus, the finite difference quotient 
of the velocity, 5u/5x, diverges with decreasing Sx, pro- 
vided Sx is still bigger than the viscous cutoff scale. Accord- 
ing to Boldyrev & Cattaneo (2004), the critical magnetic 
Reynolds number increases with increasing roughness. 

In all situations that have been simulated, the wavenum- 
ber range of the spectra has been too limited so that they 
are affected by cutoff effects both at large and small scales. 
In particular, only in simulations beyond 1024 3 meshpoints 



the spectra are shallower than fc 5 / 3 . This is referred to as 
the bottleneck effect and is believed to be a physical effect 
(Falkovich 1994, Dobler et al. 2003, Frisch et al. 2008). One 
reason, however, why it is not usually seen in wind tunnel or 
atmospheric boundary layer turbulence is the fact that one 
measures in these cases only one-dimensional spectra. In 
order to obtain three-dimensional spectra, one has to differ- 
entiate those data, i.e. (Dobler et al. 2003) 



-E-3D = -cLEm/dlnfe. 



(1) 



Accepting thus the physical reality of the bottleneck effect, 
it becomes plausible that the critical magnetic Reynolds 
number for the onset of small-scale dynamo action reaches 
a maximum around PrM = 0.1, and that it decreases some- 
what for smaller values of PrM - This is indeed what the sim- 
ulations of Iskakov et al. (2007) suggest. 

Let us now switch to large-scale dynamos. Their exci- 
tation conditions are characterized by the dynamo number 
which, for helical turbulence and in the absence of shear, is 
just 



C a — 



a 



fcf 

fcl 



(2) 



Here, fci = 2tt/ 'Lis the minimal wavenumber in the domain 
of size L and we have inserted standard approximations for 
the a effect, a = \tw ■ u, and the turbulent magnetic dif- 
fusivity, % = \tu 2 . Here, u — U — U is the fluctuating ve- 
locity, i.e. the difference between the actual velocity U and 
the mean velocity U, t ~ (u rms fcf) _1 is the turnover time, 
w = V x u is the fluctuating vorticity, ef = w ■ u/k{U 2 
is a measure for the relative helicity, and i = 1 + 3 /ReM 
is a correction factor of order unity for sufficiently large 
values of ReM- It turns out that in all cases the spectra of 
magnetic energy are at the largest scale approximately in- 
dependent of ReM for PrM between 1 and 10~ 3 . This was 
shown in Brandenburg (2009) and will here be extended to 
10 < Pr M < 10 3 . 

At larger wavenumbers there is a striking difference in 
the magnetic energy spectra between PrM = 1 and <C 1 in 
that the resistive cutoff wavenumber moves toward smaller 
values. At the same time, the kinetic energy spectrum be- 
comes progressively steeper, leaving less kinetic energy to 
dissipate. This has two important consequences. First of all, 
the fractional kinetic energy dissipation decreases with de- 

1 12 

creasing PrM proportional to Pr M (Brandenburg 2009). On 
the other hand, the decrease of 6k implies that the demand 
for numerical resolution becomes less stringent. This, in 
turn, means that one can increase the value of Re beyond the 
normally established empirical limits. An important goal of 
the present paper is the demonstration that the same is also 
true in the opposite limit of PrM 3> 1- 

3 The model 

Our model is similar to that presented in Brandenburg 
(2001, 2009), where we solve the hydromagnetic equations 
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for velocity U, logarithmic density In p, and magnetic vec- 
tor potential A for an isothermal gas in the presence of an 
externally imposed helical forcing function /, 
dU 



at 

d In p 



-U-VU-4V\xxp+f+{JxB+V-2pvS)/p, (3) 



= -U Vmp- V-U, 
U x B porjj. 



(4) 



(5) 



at 

dA 
~dt 

Here, B = V x A is the magnetic field, J = V x B/p 
is the current density, po is the vacuum permeability, c s is 
the isothermal speed of sound, and = §(J/-y + — 
■ U is the traceless rate of strain tensor. We con- 
sider a triply periodic domain of size L 3 , so the small- 
est wavenumber in the domain is k\ = 2n/L. The forc- 
ing function consists of eigenfunctions of the curl operator 
with positive eigenvalues and is therefore fully helical with 
/ • V x / = kf 2 , where 3.5 < kjk\ < 4.5 is the chosen 
wavenumber interval of the forcing function, whose aver- 
age value is referred to as fcf ss 4 k±. The amplitude of / 
is such that the Mach number is u rms /c s « 0.1, so com- 
pressive effects are negligible (Dobler et al. 2003). As in 
Brandenburg (2009), we choose as initial conditions a Bel- 
trami field of low amplitude. The initial velocity is zero and 
the initial density is uniform with p = po = const, so the 
volume-averaged density remains constant, i.e., (p) = po. 

In our simulations we change the values of magnetic and 
fluid Reynolds numbers, 

Re M = u ims /r]k{, Re = u rms /^fc f , (6) 

such that the ratio ReM /Re = Pim has the desired value 
between 10~ 3 and 10 3 , and we monitor the resulting kinetic 
and magnetic energy dissipation rates per unit volume, 

e K = (2vpS 2 ), e M = (?7MoJ 2 ), (7) 
whose sum, et = (k + £m, will be used to define the frac- 
tional dissipation rates, Ik — ^k/^t and £m — ^m/^t- 
We use the fully compressible PENCIL CODE for all our 
calculations. We recall that, for the periodic boundary con- 
ditions under consideration, (2S 2 ) = (W 2 ) + ±((V -U) 2 ), 
highlighting thus the analogy between W = V x U and J 
in the incompressible case. 




Fig. 1 Visualization of U y and B y on the periphery of the 
computational domain for Pim ranging from 10 to 1000 at 
a resolution of 256 3 mesh points. 



We note that in all cases the total energy dissipation is 
approximately the same. This is perhaps not so surprising, 
because we keep the amplitude of the forcing function the 
same. However, the constancy of the energy dissipation rate 
implies that the rate of energy injection must also be always 
the same and thus independent of PrM- This means that the 
flow properties of the eddies at the energy-carrying scale 
must be essentially independent of PrM- 

In Fig. Q] we present visualizations of the y component 
of velocity and magnetic field at the periphery of the com- 



4 Results 

In Table Q] we summarize the parameters of runs with PrM 
between 10 ~ 3 and 10 3 . The runs with 10~ 3 < PrM < 1 are 
those presented already in Brandenburg (2009) using 512 3 
mesh points, while those with 10 < PrM < 1000 are new 
ones and have been performed using 256 3 mesh points. In 
all cases, either Re or ReM were close to the maximum pos- 
sible limit at a given resolution. Indeed, for PrM = 10~ 3 we 
were able to reach Re = 4400 (for 512 3 mesh points) while 
for Pr M = 10 3 we could go to Re M = 1200 (for 256 3 mesh 
points). 

1 http ://w w w.pencil-code .googlecode.com 



Table 1 Summary of import input and output parameters 
for the runs reported in this paper. 
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Fig. 2 Compensated kinetic and magnetic energy spectra 
in the saturated regime for Pim = 10~ 3 to 10 3 . The spec- 
tra are compensated by e^ 2 ^ 3 fc 5 / 3 , where ex is the sum of 
kinetic and magnetic energy dissipation rates. The ohmic 
dissipation wavenumber, k v — (ejvf /fy 3 ) 1 / 4 , is indicated by 
an arrow. 



putational domain for the new results with Pr^ > 10 and 
in Fig. [2] we show spectra of kinetic and magnetic ener- 
gies, E(k) and M(k), respectively, for all values of PrM 
between 10~ 3 and 10 3 . In the velocity pattern one can 
clearly make out the typical scale of the dominant eddies, 
whose wave length is about 1/4 of the size of the box. 



The magnetic field also shows a turbulent component, but 
there is a much stronger large-scale component superposed. 
This is essentially the Beltrami field which is of the form 
B = (cos kiz, sinfciz, 0), although its wavevector could 
have pointed in any of the other two coordinate directions, 
(0, cos kix, sin k\x) and (sin k\y, 0, cos k\y) would have 
been equally probably alternatives. We recall that all these 
fields are indeed the eigenfunctions of an a 2 dynamo prob- 
lem (e.g., Brandenburg & Subramanian 2005), and they also 
emerge as the dominant field in helically driven turbulence. 
It is clear that in a triply periodic domain such as that con- 
sidered here, these fields require a resistive time to reach 
full saturation. For all further details we refer to Branden- 
burg (2001), where such a system was studied in full detail. 

Next, we consider the spectra of kinetic and mag- 
netic energies in Fig. [2] which are normalized such that 
J E(k)dk = \(pU 2 ) and J M(k)dk = \{B 2 /^). It is 
evident from the spectra that with increasing values of Pr^ , 
the viscous dissipation wavenumber, k v — (e/f/^ 3 ) 1 / 4 , 
moves to smaller and smaller values. Analogously to the 
case of Ptm <C 1, this implies that most of the injected 
energy gets dissipated by the shorter of the two cascades 
- leaving only a reduced amount of energy for the other 
cascade. This means that the corresponding diffusion coef- 
ficient can be decreased further, without creating numerical 
difficulties. 

It appears that it is not only the energy input at the small 
wavenumber end of the relevant cascade that is decreased, 
but that there is possibly a continuous removal of energy 
along the cascade, making the spectral index slightly steeper 
than —5/3. For example, for PrM = 10 -3 the spectral slope 
of E{k) is about —2.2, while for Ptm = 10 3 the spectral 
slope of M(k) is about -2.0. 

It is quite extraordinary that in all these cases the na- 
ture of the large-scale dynamo is virtually unchanged, even 
though Prjvr is varied by 6 orders of magnitude. The rea- 
son is that in all cases the dynamo number, C a , exceeds 
the critical value for dynamo action, C™ 1 * = 1. Looking at 
Eq. (O, we see that C a is dominated by the scale separation 
ratio, which is here kf/ki rs 4. Furthermore, because the 
turbulence is nearly fully helical, we have e{ 1, and since 
ReM 1, we have u rj 1. Thus, we have C a > 1 for all 
runs. We recall also that the saturation amplitude of the field 
is essential given by the square root of the scale separation 
ratio (Brandenburg 2001), which is about 2 in units of the 
equipartition field strength. This is in reasonable agreement 
with the simulation results; see Fig. [2] where we show the 
resulting spectra for all the runs. 

Next, we plot in Fig. [3] the ratio of kinetic to mag- 
netic energy dissipation rates. In agreement with Branden- 
burg (2009), we find that the ratio is approximately pro- 

1 /2 

portional to Pr A ^ , although a better fit is now provided by 
ejf/fM w 0.6 Prj^ 6 . The reason for such a scaling is un- 
clear. However, from Eq. (0 one can see that in the ratio 
£k /cm there is an implicit proportionality with respect to 
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Fig. 3 Dependence of the ratio of the dissipation rates on 
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where we have assumed that n lies between 1/2 and 2/3, 
which bracket the results seen here and in Brandenburg 
(2009). These scalings are surprising in view of the usually 
expected individual scalings, namely W rras oc v~ x l 2 and 
</rms °c ?7 _1 / 2 (cf. Brandenburg & Subramanian 2005). 

In order to illuminate the issue further, we ask whether 
not only the ratio ex /^m scales with PrM, but whether ex 
and em are individually proportional to Re and ReM, re- 
spectively. In Fig. 2] we plot ex versus Re (blue, solid sym- 
bols) and £m versus ReM (red, open symbols). The scatter 
is now much larger than in Fig. [3] and it seems that the scal- 
ing exponent might even be as large as n = 2 /3. 

We mentioned earlier that the total dissipation rate, ex, 
is nearly independent of PrM- However, this is only true 
when we look the the dimensional value of e-r- It is custom- 
ary to consider the normalized dissipation rate, 

c t = tT 

u 3 1B /L 

where mid = u rms /^/3 is the one-dimensional rms velocity 
and L = 37r/4fcf is conventionally used as the integral scale 
(Pearson et al. 2004). In the second and third panels of Fig. [4] 
we compare C e with C e o, which is based on the maximum 
value of uid in all the runs. The difference is caused by the 
fact that u Tms drops to rather low values in the large-PrM 
regime. Part of this goes into magnetic energy, but it is not 
enough to make up for this difference. 

It is important to realize that, on average, eu is just the 
same as the rate of work done against the Lorentz force, 
— {U ■ (J x B)). This becomes evident when considering 
the flow of energy in our system: 

-> (2puS 2 ) 
-(U-(JxB)) ^( W oJ 2 
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Fig. 4 Top: Dependence of ex on Re (blue, solid sym- 
bols) and £m on ReM (red, open symbols). The solid line 
has the slope 2/3, while the dotted and dashed lines have 
slopes 0.6 and 0.5, respectively. Middle and bottom: scal- 
ings of C e o and C t versus PrM- 



Here, (pU ■ f) ~ e t is the rate of energy injection into the 
system by the forcing term. Normally, in the hydrodynamic 
( ' case, (2pvS 2 ), or (vW 2 ) in the incompressible case, stay 
constant as v is decreased. In the case with dynamo action, 
however, a decrease in v allows the dynamo to tap more 
energy, so — (U ■ (J x B)) and cm increase at the expense 
of ex. This is indicated by the fact ex /^m is found to be 
proportional to (y/rj) n , so ex decreases as v decreases. This 
decease is weak in the sense that n « 1/2 ... 2/3 is less than 
unity, but it is certainly no longer independent of v as it 
would be in the purely hydrodynamic case. 

In view of the application to quasars, i.e. accretion discs 
in active galactic nuclei, it is relevant to consider the frac- 
tion of energy that goes into the heating of electrons. In- 
deed, such discs are known to be underluminous, which led 
to the standard paradigm of advection-dominated accretion 
^ ' (Narayan & Yi 1994; Abramowicz et al. 1995). Alterna- 
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tively, this might be associated with the small value of the 
ratio em /et, for which we find 



1 



oc 



e K 1 + M 



(12) 



Using standard accretion disc theory, Balbus & Henri 
(2008) find that Pim depends on the distance R from the 
black hole and is proportional to i?~ 9 / 8 . In particular, they 
find that PrM exceeds unity within about 50 Schwarzschild 
radii. This would dramatically decrease cm in the inner 
parts and might be sufficient to explain underluminous ac- 
cretion. However, this proposal hinges on several assump- 
tions: (i) that the viscous heating heats the ions and not the 
electrons, (ii) that the resistive dissipation energizes elec- 
trons rather than ions, (iii) that the discs are essentially col- 
lisionless and, finally, (iv) that the magnetohydrodynamic 
approximation is then still applicable. 



5 Conclusions 

The present work has shown that the ratio of kinetic to mag- 
netic energy dissipation follows one and the same relation- 
ship with PrM both for small and large values. An impor- 
tant additional condition obeyed by all our runs is, however, 
that the magnetic Reynolds number is large enough for dy- 
namo action to occur. This constitutes an important differ- 
ence between our current results for large-scale dynamos 
and those mentioned in the first section for small-scale dy- 
namos. An important consequence of such scaling is the 
fact that at extreme values of the magnetic Prandtl number, 
larger Reynolds numbers can be tolerated by the numerical 
scheme at a resolution that would be insufficient if the mag- 
netic Prandtl number were unity. This was shown previously 
for PrM = 10 -3 , in which case fluid Reynolds numbers 
of 4500 were possible at a resolution of 512 3 meshpoints, 
while for PrM = 1 it was only possible to reach Reynolds 
numbers of less that 700. Both cases obeyed the empirical 
constraint that the spectral kinetic energy has developed a 
clear dissipative subrange with an exponential decay shortly 
before the Nyquist frequency. In the opposite case of large 
PrM, here PrM = 10 3 , it was possible to reach magnetic 
Reynolds numbers of 1000 at 256 3 mesh points. In this case 
the magnetic energy spectrum has developed a dissipative 
subrange shortly before the Nyquist frequency, although it 
was less convincing for PrM = 10 2 . 

The reason for the value of the exponent n in the power 
law relation between the energy dissipation ratio ck/^m 
and PrM remains unclear. At this point we cannot be cer- 
tain that it is n = 0.6 and not, for example, 1/2 or 2/3. 
One source of error might come from the fact that at ex- 
treme values of PrM the effects of numerical viscosity as- 
sociated with the advection operator are no longer negligi- 
ble. For the third-order time step used in the PENCIL Code, 
the numerical viscosity operator takes the form — z/g L V 4 
6x 3 Cq FL /24 is a numerical hypervis- 



cosity that depends on the mesh size Sx and the Courant- 
Friedrich-Levy number Ccfl, whose default value is 0.4, 
but the code would still be numerically stable for Ccfl = 
0.9. If such numerical effects do begin to play a role, we 
must expect that the effective values of PrM are less ex- 
treme, which means that the n would have been underesti- 
mated and that n might be 2/3 or even larger. 

While earlier work focussed on the dependence of £m 
on PrM (Blackman & Field 2008), no clear conclusion 
about the dissipation ratio ck Am seems to have emerged. 
For example, if ejf and cm were independent of viscosity 
and magnetic diffusivity, the ratio ck /^m would have been 
constant. Instead, we find that ex decreases when Re de- 
creases, and likewise, cm decreases when ReM decreases. 
On the other hand, one must be cautious when applying re- 
sults regarding the dependence on Rejif/Re (= Pr^) for 
large values of Re and ReM, because we may still not be in 
an asymptotic parameter regime. It is therefore important to 
extend this work to larger values of Re and ReM and to go 
to larger numerical resolution. 
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where ^2 !FL = u 



1 See page 118 of the Pencil Code manual, 

http : // www . nordita . org/ software/pencil- code /doc /manual .pdf 
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